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We consider semi-infinite two-dimensional layered Ising models in the extreme anisotropic limit 
with an aperiodic modulation of the couplings. Using substitution rules to generate the aperiodic 
sequences, we derive functional equations for the surface magnetization. These equations are solved 
by iteration and the critical exponent (is can be determined exactly. The method is applied to 
three specific aperiodic sequences, which represent different types of perturbation, according to a 
relevance-irrelevance criterion. On the Thue-Morse lattice, for which the modulation is an irrel- 
evant perturbation, the surface magnetization vanishes with a square root singularity, like in the 
homogeneous lattice. For the period-doubling sequence, the perturbation is marginal and Ps is a 
continuous function of the modulation amplitude. Finally, the Rudin-Shapiro sequence, which cor- 
responds to the relevant case, displays an anomalous surface critical behavior which is analyzed via 
scaling considerations. Depending on the value of the modulation, the surface magnetization either 
vanishes with an essential singularity or remains finite at the bulk critical point, i.e., the surface 
phase transition is of first order. 



I. INTRODUCTION 

Since the discovery of quasicrystals,S one has witnessed 
a growing interest in understanding their structure and 
physical properties (for recent reviews see Refs. |^-|^). 
The study of phase transitions on quasiperiodic lattices is 
now an active field of research. Since quasiperiodic and, 
more generally, aperiodic systems can be considered as a 
state of matter which interpolates between the periodic 
(crystalline) and the random (disordered, glassy) states, 
the corresponding critical behavior is very rich and the- 
oretically challenging. 

The first results obtained in this field were mainly 
numerical. For example, it was shown that different 
systerna jOn the two-dimensional Penrose lattice (Ising 
modeljErEl percolationJlBO self-avoiding walkst3) are uni- 
versal in the sense that the critical exponents are the 
same as in regular, two-dimensional lattices. The same 
conclusion was jEeached for three-dimensional quasiperi- 
odic structures ,lJ although with a weaker numerical ac- 
curacy. 

Exact results about critical properties can be obtained 
on two-dimensional, layered Ising models, where the in- 
tralayer couplings Ki are kept constant, whereas the in- 
terlayer interactions K2{k) follow the aperiodic modu- 
lation of the underlying lattice. One usually works in 
the extreme anisotropic limit JU with K2{k)^Q, Ki^oo, 



keeping fixed the ratio Afc — K2{k) / Kl where Kl is the 
dual interaction corresponding to Ki. In this limit the 
system is described by a one-dimensional quantum Ising 
model (QIM) with the Hamiltonian 

^ = -^Ek' + ^fc^fc^fc+i]' (1-1) 
fc=i 

where the tr's are Pauli spin 1/2 operators. 

The first exact results on the QIM were obtained for 
the bulk critical properties on the Fibonacci and related 
lattices, where the specific heat is found te-have a log- 
arithmic singularity of the Onsager typeJl3"llj It was 
Tracyc^ who first showed examples of aperiodic layered 
Ising models where the logarithmic singularity in the spe- 
cific heat is washed out, like in the layered -laing model 
with random couplings (McCoy- Wu model) 

A systematic study of the bulk critical behavioij-of ape- 
riodic QIMs has been performed-recently by LuckB3 Gen- 
eralizing the Harris criterion,Ej the relevance or irrele- 
vance of the aperiodicity is shown to be connected to the 
size of the fluctuations in the couplings Afc . For bounded 
fluctuations (as happens for the Fibonacci lattice) the 
speciflc heat of the system diverges logarithmically as 
in the homogeneous one, while, for unbounded fluctua- 
tions, the critical behavior is anomalous, e.g., the specific 
heat displmis an essential singularity like in the McCoy- 
Wu model. Eil In the marginal case, when the fluctuations 
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grow on a logarithmic scale, non-universal critical behav- 
ior is expected with critical exponents depending on the 
strength of the aperiodicity. 

The relevance- irrelevance— .criterion has been gen- 
eralized for Qihei systemsEil and higher dimensional 
aperiodicitiesBa For relevant modulations, the form of 
the singular quantities near the critical ueint has also 
been determined, using scaling arguments.cJ 

As far as the surface critical behavior of aperiodic sys- 
tems is concerned, only a few results are available. For 
the QIM in the ordered phase, the asymptotic limit of 
the spin-spin correlation function in the surface gives the 
square of the jSurface magnetization rris . This leads to 
the expressionE3 



nis =< l|o-i |0 > 



(1.2) 



where |0 > is the ground state of H and |1 > is the first 
excited state, belonging to the odd sector of the Hamilto- 
nian, which is degenerate with the ground state in the or- 
dered phase of the infinite system. The matrix element in 
Eq. (1.2) can be rewritten using a Jordan- Wigner trans- 
formation of the spin operator, followed by a canonical 
transformation to Fermi operators which diagonalizes the 
Hamiltonian as described in Ref. |2^. The surface mag- 
netization finally takes ai-simple form involving sums of 
products of the couplinggHl 



'1/2 



rus 



3 = 1 k=l 



(1.3) 



which stays valid for any distribution. The surface mag- 
netization has been evaluated in the critical region for 
different sequences with bounded fluctuations in Ref. |2^ 
and in each case it was found to vanish with a square law 
singularity, in agreement with scalin g a rguments. On 
the other hand it was shown in Ref. ^ that any type 
of marginal modulation results in a non-universal critical 
behavior, with a surface magnetization exponent which 
is a continuous function of the modulation strength. 

In the present paper, we continue and extend our stud- 
ies of the surface critical behavior of aperiodic quantum 
Ising models. To evaluate the infinite sum for in 
Eq. (L3) we have developed a method leading to fimc- 
tional equations for the surface magnetization which do 
not contain explicitly the form of the aperiodicity. It- 
erating these relations, a closed form expression for the 
surface magnetization is obtained, from which the critical 
exponent as well as corrections to scaling can be deter- 
mined exactly. To illustrate the method, different aperi- 
odic sequences with irrelevant, marginal, as well as rele- 
vant modulations, are evaluated. For the Rudin-Shapiro 
sequence, which represents a relevant perturbation, a 
combination of numerical results and scaling arguments 
has been used. The surface magnetization is found to 
display a first order transition for some range of the cou- 
pling ratio, whereas the bulk transition is expected to be 
continuous. 



The setup of the paper is the following. In Sec. ||, 
the properties of aperiodic sequences generated through 
substitutions are summarized and a relevance-irrelevance 
criterion is deduced from scaling considerations. Then we 
study successively the Thue-Morse sequence (Sec. p| ), 
the period-doubling sequence (Sec. IV), and the Rudin- 
Shapiro sequence (Sec. The results are discussed in 
the final section. 



II. APERIODIC SEQUENCES AND SCALING 
CONSIDERATIONS 

Aperiodic sequences can be generated through iterated 
substitutions on the letters A, B, ... such that A 
S{A), B~iS{B), .... The properties of a given sequence 
are then obtained from its substitution matrix M whose 
columns contain the numbers of letters A, B, . . . in S{A), 
S{B), . . . , respectively: 



M = 



S{A) 5(B) 
S{A) SiB) 



'B 



'B 



(2.1) 



It follows that M" has matrix elements giving the corre- 
sponding numbers in the sequences constructed on A, 
B, ... after n substitutions. For example, when one 
starts with A, i.e., with the substitutions A S{A) 
S {S{A)) . . ., the number of A, L^, in the sequence after 
n steps is (M"')^^ whereas the length of the sequence is 
Ln — J2i {MDii ■ Making use of the right eigenvectors Va 
of the substitution matrix 



(2.2) 



Ln and are found to be asymptotically proportional 
to A", where Ai is the largest eigenvalue of the substitu- 
tion matrix, whereas the asymptotic density of A letters 
involves the components of the corresponding eigenvec- 
tors 



Via) 



Fluctuations in the numbers of letters L„, L„ , 



(2.3) 



are 



connected to the next-to-leading eigenvalue A2. 

In the case of the aperiodic quantum Ising chain con- 
sidered above, with a two-letter sequence {Xk = Aa, As) 
and an averaged coupling A, one may write 



Xa 



Xb=X 



A 1 



(2.4) 



where S ~ Xa — Xb- Then, at a length scale Ln ^ A^, 
the sum of the deviations from the averaged coupling is 



A(L„) 



E 

k=l 



{Xk - A) 



SL't 



(2.5) 
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and involves the wandering exponent"''^'^'P^' 
ln|A2| 



(2.6) 



When the number of letters is greater than 2, one may 
have two complex conjugate next-to-leading eigenvalues 
and the power law in Eq. ( ^.5| ) is modified by a factor 
which is periodic in n ~ lnL„/lnAi. 

The aperiodicity introduces a thermal perturbation 
above A which, at a length scale L, is given on the average 
by 



(2.7) 



Under a change of scale by a factor h = L/L' this trans- 
forms according to 

W{L') = S' (^^^ = b^^''SL^-\ (2.8) 

where v is the bulk correlation length exponent. Jtjpl- 
lows that the perturbation amplitude S scales likeESEJ 



S' = h^l^i, 
with a crossover exponent 

$ = 1 + V{'jj - 1). 



(2.9) 



(2.10) 



As a consequence, for the two-dimensional Ising model 
with v — \, the aperiodic modulation is a relevant (irrel- 
evant) perturbation when cj > (tj < 0) and it becomes 
marginal when lo = [}. In the latter case a nonuniversal 
behavior is expected. 

For a relevant aperiodic modulation, the form of the 
singular thermodynamical quantities can be deduced 
from scaling considerations. Let us consider the surface 
magnetization as a function of the thermal scaling field 
t = 1 — (Ac/A)^ and the modulation amplitude b. In a 
scale transformation it behaves as 

m,(t, S) = h-^^l^'m, (b^/^t, , (2.11) 

where Ps is the surface magnetization exponent. Taking 
the scale factor b equal to the bulk correlation length 
b = ^ = t~'^ , one obtains in the form 



(2.12) 



where the scaling function F{x) involves a new charac- 
teristic length. 



(2.13) 



which is introduced by the aperiodicity and remains finite 
at the bulk critical point. 



In the following, an aperiodic sequence will be writ- 
ten as a succession of digits fk = 0, 1, which are im- 
ages 3(A, B) of the letters considered before. Then, with 
Ai — Ar, Ao = A, one obtains A^ = Ar^*" for the fcth cou- 



pling and the surface magnetization in Eq. (1.3) can be 
rewritten as 



where 



3 

k=l 



5(A,r)=^A- 



no = 0. 



(2.14) 



(2.15) 



The critical coupling, A = A^,, is such thatE 
limL^oo rifc^i Afcc = hmL^oo A^r"^ =1, so that 



= lim — ^ . 

L — ^oo L 



(2.16) 



III. THUE-MORSE SEQUENCE 

As a first example, we consider the binary Thue-Morse 
sequence,L3 which may be deduced from the binary rep- 
resentation of non-negative integers, 0, 1, 10, 11, . . ., by 
counting the sum of their digits modulo 2, 

01101001 ... . (3.1) 

This sequence is also generated through the substitution 



AB 
BA 



1 

1 



(3.2) 



with 9(A, B) — (0, 1), so that one obtains, successively, 


1 

110 

110 10 1 

0110100110010110 

(3.3) 

The corresponding substitution matrix ( p.l[ ) has eigenval- 
ues Ai = 2, A2 = 0. The asymptotic density, Po o — 1 /2, 
is deduced from Eq. (2.3) and, according to Eq. (2.16|), 



A. 



-1/2 



(3.4) 



The form of the substitution in Eq. (3.2) immediately 
leads to the recursion relations 



/2p - 1 - /i 



Pi 



/2p+l — /p+1- 



(3.5) 



From the sec ond one, the sequence of odd terms, under- 
lined in (3.3), reproduces the whole sequence. Splitting 
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t he sum giving rij into even and odd parts and using 
( |3.5|) , one obtains 



n2p=P, n2p+i = p+ fp+i, p = 0,1,2,..., 

(3.6) 

so that a chain with length L ~ 2p has a density equal to 
the asymptotic one, poo = 1/2; which explains the vanish- 
ing second eigenvalue giving oj — —oo, i.e., no wandering 
in this case. 



0.8 



« 0.6 



0.4 



0.2 




FIG. 1. Thue-Morse surface magnetization as a function 
of l — t— (Ac/A)^ for different values of the coupling ratio 
r = Ai /Aq. The critical exponent keeps its value for the homo- 
geneous Ising system /3s = 1/2 whereas the critical amplitude 
is r-dependent. 



Th e sum S{X, r) involved in the surface magnetization 
(2.14) can be rewritten as 



5(A,r)=^A-2(2p)r-2n2p 



-2(2p+l) -2n2p+i 



p=0 
oo 



p=0 



2(2p+l)„-2p-2/p+ 



K (3.7) 



p=0 



p=0 



Using the identity a^'' = l + {a— l)fk jfk — 0, 1) and the 
value of the critical coupling in Eq. (3.4), one obtains 



1 - (Xc/xy 



(3.8) 



fc=i 



The Thue-Morse series, S(a: ), sa tisfies a recursion rela- 
tion which follows from Eq. (|3.5|), 



^x) 



1 — X' 



+ {x-^ - l)S(a;2), 



(3.9) 



so that, iterating this relation. 



S(a;) 



^ (l-a;2M (1 
fe=o ^ > ^ 



^2fe+n 



(3.10) 



whic h is c onvergent for a; < 1. Equations ( |2.14 ), ([3.^)^ 
and ( |3lc| ) give the surface magnetization shown in Fig. nJ 



for any value of A > Ac 
Near the critical point, A 



Ac+ , the leading contribu- 



tion to ^(A, r) comes from the first two terms in Eq. (p.q). 



The next ones in E 



(Ac/A)^ 



give the corrections to scal- 



ing. Collecting these results, the surface magnetization 
finally takes the form 



Ac + Ac 



I -XI 
1 + A2 



t + 0[t^ 



(3.11) 



where t is the deviation from the critical point defined 
in Sec. The surface magnetization exponent takes 
the value Ps = 1/2 for homogeneous Ising systems in 
two dimensions as expected from scaling since u; < 
here and the aperiodic modulation of the couplings is an 
irrelevant perturbation. The critical amplitude depeads 
on r through Ac in agreement with previous results.E3 



IV. PERIOD-DOUBLING SEQUENCE 



We now turn to the pericjdj-doubling sequence which 
follows from the substitutionEJO 



1 
1 1 



(4.1) 



leading to 



1 

1 

10 11 

10 1110 10 

1011101010111011 



(4.2) 



The eigenvalues of the substitution matrix are then Ai = 
2, A2 = — 1 and the asymptotic density poo = 2/3 gives 
the critical coupling 



Ac = r-2/3. 



(4.3) 



The form of the substitution is such that /2p = 1 — ./p, 
/2p+i = 1, from which one deduces 



n2p = 2p - Up, n2p+i ^2p+l 



p = 0,l,2, 



(4.4) 



Splittin g th e sum in S( X, r) into even and odd parts as 
in Eq. (3/7) and using {iA), one obtains the following 
recursion relation: 
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5(A,r)= (l + -f^)5(AV,r-i) 



(4.5) 



After I iterations the arguments in S{X, r) are changed 
into 



A2Vi(2'+i), / odd, 



r^'^ ^\ I even, 

^ = 0,1,2,..., (4.6) 

and the series can be written as an infinite product 



(-1) 
ri— ' 



5(A,r) = n 



fe=i 



(4.7) 



which, together with (|2.14 ), gives the surface magnetiza- 
tion shown in Fig. ^ 




(kA) 

FIG. 2. Period-doubling surface magnetization for different 
values of r. The exponent l3s varies continuously with the 
modulation amplitude. 

In order to analyze the critical behavior of the sur- 
face magnetization, jjMie may use a scaling method intro- 
duced by one of usO Let S{u) denote the series expan- 
sion of 5' (A, r) in powers of m = (Ac/A)^. According to 
Eq. ( 2.14 ), S{u) should display a power law singularity 
near the critical point: 



S{u) - {l~u) 



-20s 



(4.8) 



It may be shown that in this case, the truncated series 
Sl{u) given by the first L terms in S{u) behaves, for large 
values of L, as L^^' at the critical point, u = 1. Now 
it may be verified that the first I terms in the infinite 
product (4.7) just contain the L = 2^' first terms of the 



series expansion. As a consequence, one obtains 

^L=2- (1) = [(1 + Ac)(l + A-i)] ' ^ (2^f 
from which the surface magnetization exponent 



(4.9) 
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4 In 2 



(4.10) 



follows. 

The same result can be obtained by considering the 
recursion relation (4.5) at the next step: 



S{\r)= 1 + 



1 



A2r2 



1 



1 



A4r2 



5(AV2,r). (4.11) 



At this stage the new coupling on the right-hand side 
is A' = A^r^ whereas r and, as a consequence, Ac re main 
unchanged. When A ^ Ac, S{\, r) behaves as in Eq. ( |4.8| ) 
with an amplitude A{r) and, in the same way. 



5(AV2,r) ~ A(r) 



A{r) 



1- ^ 



-2/3, 



-2/3, 



(4.12) 



where the value of Ac given in Eq. (4.2) was used. In- 
troducing these results i nto ( 4.11) an equation for [3s is 
obtained which leads to (4.1Cj! 

The surface magnetization exponent depends on the 
amplitude of the aperiodicity r through Ac as expected 
for this sequence since u vanishes and the perturbation 
is marginal. The variation is shown in Fig. ^. 



CO. 
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FIG. 3. Marginal variation of surface exponent /3s with the 
coupling ratio r for the period-doubling sequence. In the ape- 
riodic system the singularity is always weaker than in the ho- 
mogeneous one and the surface transition is of second order 
for any value of r. 



Let ms{L) be the matrix element in Eq. (1.2) for a 
chain with length L. Although the spontaneous magne- 
tization vanishes on a finite chain, ms{L) remains ncpyja-. 
nishing and scales with L as L~^= at the critical pointEa^ES 
where Xs = (3s I v is the scaling dimension of the surface 
magnetization . To the leading order in L, ms{L) is still 
given by Eq. (L3) with a sum over j — I, L so that the 
truncated series Sl{1) considered above also behaves as 



L. TURBAN, F. IGLOI AND B. BERCHE 



Phys. Rev. B 49 (1994) 12700 



L^^". It follows that Ps = Xs and v =1: The marginal 
aperiodic modulation does not change the bulk correla- 
tion length exponent. 



V. RUDIN-SHAPIRO SEQUENCE 



As a 
sequence,! 
stitution 



example, we consider the Rudin-Shapiro 
i3y22| jg generated using a four-letter sub- 



< 



A 
B 
C 



AB 
AC 
DB 
DC 



(5.1) 



One may take eithat, single-digit images of the four 
letters in (|5T|) withB ^A,B,C,D) = (0,0,1,1) or, 
more naturally, two-digit images with Q[A, B,C, D) = 
(00,01,10,11) so that 



00 ^ 0001 

01 ^ 0010 
10 1101 
11 1110 



(5.2) 



The two processes lead to the same sequence of 1 and 
since the two-digit images are related to the single-digit 
ones through a substitution on the letters. Starting on 
A, one obtains 



1 

10 10 

0001001000011101 

(5.3) 

The substitution matrix has leading eigenvalues Ai = 2, 
A2 = ±\/2 so that w = 1/2 and the asymptotic density 
Poo = 1/2, obtained by taking a weighted average on 
the densities for the different letters, gives the critical 
coupling 



A. 



-1/2 



(5.4) 



From the form of the substitution in Eq. (5.2), one easily 
deduces the following relations: 



/4p — 1 — f2p, /4p+l — /4p+2 — /: 



.fip+3 — /2p+2- 

(5.5) 



Furthermore, one may notice that odd terms, underlined 
in (^) , reproduce the sequence itself whereas even terms 
either reproduce the sequence or its complement to one, 
so that: 



2p+l 



fp 



+ I1 



f: 



2p+2 



1 - /p+1, 



p odd, 
p even. 



(5.6) 




(kA) 



FIG. 4. Surface magnetization for the Rudin-Shapiro se- 
quence with a coupling ratio r < 1. The enhancement of the 
coupling strength near to the surface leads to a first-order 
surface transition when r < 1. 



« 0.6 




FIG. 5. Surface magnetization for the Rudin-Shapiro se- 
quence with r > 1. The surface magnetization vanishes with 
an essential singularity at the bulk critical point when r > 1. 

Using these p rope rties one may write a recursion rela- 
tion similar to (4^), but here in matrix form, between 
the odd and even parts of S{X,r). Although the matrix 
considerably simplifies at the bulk critical point, some of 
its elements still depend on the iteration index /. As a 
consequence, we were unable to deduce the critical be- 
haviour from this recursion relation like in Sec. 



[V 



Numerical results for the surface magnetization are 
shown in Figs. ^ and |5| for r < 1 and r > 1, respec- 
tively. In the first case, one obtains a first order surface 
transition, which is linked to a local enhancement of the 
couplings relative to their averaged value. On the other 
side, when r > 1, the couplings near the surface are lo- 
cally weaker than the average and the surface magnetiza- 
tion vanishes at the bulk critical point with a singularity 
which is seemingly weaker than any power law. In the 
bulk, the situation is different since the deviations from 
the average have both signs with the same probability. 
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That the first order transition is indeed a surface effect 
can also be made plausible through the following argu- 
ment: If one started the sequence with D instead of A, 
the bulk properties should not change but one can check 
on Eqs. (5.1) and (5.2) that it would amount to exchange 



weak and strong couplings along the sequence, and the 
first order surface transition would then occur for r > 1. 

The discontinuity of the critical surface magnetization 
shown in Fig. ^ appears to vary linearly with the coupling 
ratio r < 1. This behavior will be explained in Sec. 
using scaling arguments. 



0.12 



VI 




1.02 



FIG. 6. Surface magnetization discontinuity at the bulk 
critical point as a function of the coupling ratio for the 
Rudin-Shapiro sequence with increasing sizes, L = 2^^^ to 
2^*, from top to bottom. The variation is asymptotically 
linear in r. 



VI. DISCUSSION 

The QIM on aperiodic lattices studied in this paper 
shows a rich critical behavior. Depending on the value of 
the wandering exponent w, scaling arguments show that 
the perturbation introduced by the aperiodic modulation 
can be either irrelevant, marginal or relevant. This has 
been verified on three specific sequences for which the 
different typical behaviors are observed. 

For the period-doubling sequence, displaying a 
marginal modulation, the surface critical exponent [3^ 
and the corresponding scaling dimension Xg = Ps/i^ are 
nonuniversal and show the same continuous variation 
with the strength of the modulation. As a consequence 
y itself stays constant, keeping its value in the homoge- 
neous system, v = 1. Thus the marginality condition 
1/u = 1 — w is fulfilled along the critical line, as required, 
On the other hand one knows from numerical studiesEHI 
that the bulk specific heat exponent a does vary with 5 
so that the hyperscaling relation dv — 2 — a is violated for 
this marginal sequence. Such a violation of hyperscaling 
is likely to occur for other marginal sequences too. It may 
be due to the pptt^ence of a dangerous irrelevant variable 
in the problem.tZI 
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FIG. 7. Logarithm of the reduced surface magnetization 
as a function of for the Rudin-Shapiro sequence with 
r > 1 and L = 2^^ . The linear variation (dashed line) is i n 
agreement with the stretched exponential form in Eq. (6.3). 
Deviations from scaling are mainly due to finite-size effects. 



With the Rudin-Shapiro sequence, the relevant modu- 
lation was found to lead to different surface critical be- 
haviors, depending on whether the coulings near to the 
surface are greater or smaller than their average. In the 
first case, the surface stays ordered at the bulk critical 
point. The variation with r of surface magnetization 
discontinuity can be deduced from the scaling form in 
Eq. ( ^.12 ). For ms(0, 5) to remain finite at the bulk crit- 
ical point, the t dependence on the right-hand side has to 
cancel. The scaling function then behaves as a power of 
its argument which follows from the cancellation, leading 
to 



m,(0,(5)^/. 



1^1 



A/* 



(6.1) 



In the present case the exponent is equal to 1 and the 
discontinuity varies as |^l = 1 — r in agreement with the 
numerical results in Fig. @. 

The same behavior is, obained in the the Hilhorst-van 
Leeuwen (HvL) modeCa in the case of an inhomogeneity 
decaying as AL~y where L is the distance from the free 
surface. In this cas e the crossover exponent is 1 — vy so 
that, according to (2.10), y corresponds to 1 — w which, 
in the present problem, is the decay exponent for the av- 
erage deviation of the couplings given in Eq. (2.7). The 
jump of the surface magnetizatian in the Ising HvL model 
has been obtained analyticallyEa for a relevant perturba- 
tion with A> Q and is in agreement with ( |6.l|) . 

When r > 1, the local couplings at the surface are in 
average weaker than the critical one and, according to the 
numerical results shown in Fig. |^, the surface order van- 
ishes in an anomalous way. In order to obtain the form 
of the scaling function F{x) in Eq. ( |2.12D , one may use 
the abovementioned analogy with the Hvlj. model as done 
in Ref. ^ An argument of Burkhardt,Ej going beyond 
scaling theory, leads to a stretched exponential behavior 
for the critical correlation function, from which the lead- 
ing singularity of the susceptibility can be deduced,L3 
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X ~ exp 



(6.2) 



In this expression x = la/S, is the scahng variable defined 
in Eq. ( p.l2| ). The same form is likely to occm- in the 
temperature dependence of other singular quantities and 
one expects the reduced surface magnetization to behave 
as 



rris 



exp 



(6.3) 



where Cm is some constant. The scaling function may 
also involve some unknown power of x in front of the 
exponential. This expression, when properly translated, 
is in agreement with the analytical [Result obtained by 
Peschel for the surface inhomogeneity.Ea For the aperiodic 
Ising model with Rudin-Shapiro modulation, lu — 1/2 
and ly — 1, so that the argument of the exponential in 
Eq. (|6.3|) is a:"'^ = (r — l)^/i. The numerical results shown 
in Fig. |7|, although perturbed by finite-size effects, con- 
firm the proposed scaling form. 

All these systems containing a characteristic length, 
which, like la, stays finite at the bulk critical point when 
the perturbation is relevant, seem to display the same 
type of critical behaviour. Besides the two examples con- 
sidered so far, one may also mention systems limited by a 
surface with a parabolic shape, for jaihich the character- 
istic length is fixed by the geometry.Ej They also differ in 
some aspects: In the marginal HvL model, for instance, 
the surface transition is first order when the amplitude of 
the inhomogeneity is positive and strong enough, whereas 
it remains second order for the period-doubling modula- 
tion studied here, whatever the modulation amplitude. 
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